Comparative experimental study of local mixing of active and passive scalars in 

turbulent thermal convection 
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We investigate experimentally the statistical properties of active and passive scalar fields in turbu- 
lent Rayleigh-Benard convection in water, at Ra ~ lO'^". Both the local concentration of fluorescence 
dye and the local temperature are measured near the sidewall of a rectangular cell. It is found that, 
although they are advected by the same turbulent flow, the two scalars distribute differently. This 
difference is twofold, i.e. both the quantities themselves and their small-scale increments have dif- 
ferent distributions. Our results show that there is a certain buoyant scale based on time domain, 
i.e. the Bolgiano time scale is, above which buoyancy effects are significant. Above tB, temperature 
is active and is found to be more intermittent than concentration, which is passive. This suggests 
that the active scalar possesses a higher level of intermittency in turbulent thermal convection. It is 
further found that the mixing of both scalar fields are isotropic for scales larger than ts even though 
buoyancy acts on the fiuid in the vertical direction. Below tB, temperature is passive and is found 
to be more anisotropic than concentration. But this higher degree of anisotropy is attributed to the 
higher diffusivity of temperature over that of concentration. From the simultaneous measurements 
of temperature and concentration, it is shown that two scalars have similar autocorrelation functions 
and there is a strong and positive correlation between them. 

PACS numbers: 47.27.-1, 47.27.Te, 47.51. +a, 47.55.P 



I. INTRODUCTION 

In the studies of hydrodynamic turbulence, scalar field 
has been one of the main subjects besides the velocity 
field. This is partly because a scalar field is presumed 
to be easier to tackle with than a vector field, both ex- 
perimentally and theoretically. In general, there are two 
main kinds of scalar fields, classified according to the in- 
teraction between the scalar and its carrier flow: One is 
active scalar which couples dynamically to the velocity 
field and the other is passive scalar whose feedback to 
the flow is negligible. Although both are governed by 
the same advection-diffusion equation, active and pas- 
sive scalars belong to different realms of mathematics 
and physics. Due to the dynamical coupling with the 
advecting velocity field, the problem of active scalars is 
a nonlinear one. On the other hand, passive scalars obey 
a linear equation because of the absence of the feedback 
to the advecting velocity field. This allows a fully theo- 
retical treatment of the problem, which has been carried 
out recently [J, • To study the properties of both scalar 
fields, the turbulent Rayleigh-Benard convection (RBC), 
aparadigm for studying buoyancy-driven turbulent flows 
0, [1, H, S, 0| I provides an ideal platform. This is because 
in this system the temperature field itself is already an 
active scalar, as it drives the turbulent flow via buoyancy, 
and, if a fluorescent dye solution is injected into the flow, 
the concentration field of the dye behaves as a passive 
scalar. The object of the present experimental investiga- 
tion is to address the following questions: What are the 
relations, similarities, and differences between the active 
and passive scalar fields advected by the same turbulent 
flow? 

Passive scalars have been studied extensively for many 



years, and the results have provided many insights into 
the turbulence problem [J, [3, B HI ■ However, there have 
been a limited number of studies on the relations between 
the active and passive scalar fields. Celani et al. [l3| com- 
pared the behaviors of active and passive scalars in four 
different numerical systems [Tol.[TllTT^ by focusing on the 
issue of universality and scalings and their results provide 
two possible scenarios: the first is that active and pas- 
sive fields in the same flow should share the same statis- 
tics if the statistical correlations between scalar forcing 
and carrier flow are sufficiently week; whereas in the sec- 
ond scenario these two fields may behave very differently 
for systems with strong correlations between the active 
scalar input and the particle trajectories. These authors 
thus suggested that a case-by-case study is needed. For 

ll| showed that 



thermal convection, Celani et al. [1 
in 2D numerical convective turbulence the 



two scalar 

fields have the same even-order anomalous correlations, 
which are universal with respect to the choice of exter- 
nal sources, and their intermittency in such systems may 
be traced back to the existence of statistically preserved 
structures. However, these authors also argued that the 
equivalence of the statistics of active and passive scalar 
fields depend crucially on the statistics of the whole veloc- 
ity field, which may not hold in 3D convective turbulence. 
Ching et al. also offered numerical evidence in the 
context of simplified shell models for turbulent convec- 
tion to show that under generic conditions the even-order 
correlation functions of active scalars can be understood 
via the emerging theory of statistically preserved struc- 
tures of the passive scalar counterpart. For magnetohy- 
drodynamics, Celani et al. [13, [13 showed that a passive 
scalar displays a direct cascade towards the small scales 
while the active magnetic potential builds up large-scale 
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structures in an inverse cascade process. Gilbert and 
Mitra further found that active and passive fields 
have different scaling properties in the framework of a 
shell model. Ching et al. argued that such differ- 
ent scaling behaviors are due to the fact that the active 
equations possess additional conservation laws and the 
zero modes of the passive problem is not the leading fac- 
tor that dominates the structure functions of the active 
field. However, some of the above fundamental relations 
between active and passive scalar fields have not been 
investigated experimentally. 

In an earlier study, Zhou and Xia 15| investigated 
experimentally the mixing of an active scalar, i.e. the 
temperature field, in thermal turbulence. Their results 
show that the active scalar field has sharp fronts or gra- 
dients with similar statistical properties as those found in 
passive scalars, such as the saturation of structure func- 
tion exponent and log-normal distribution of the width of 
fronts. However, to the best of our knowledge, there have 
been no experimental investigation to directly compare 
the fundamental properties of active and passive scalars 
that are driven by the same turbulent velocity field. In 
this paper, we undertake such studies in a turbulent RBC 
system and our results suggest that the active scalar pos- 
sesses a higher level of intermittency. 

The remainder of this paper is organized as follows. In 
Sec. ini we describe the convection cell used in the exper- 
iments, details of the temperature measurement and the 
technique of laser-induced fiuorescence (LIF), and exper- 
imental conditions and parameters. The experimental 
results are presented and analyzed in Sec. IIIII which 
is divided into five parts. Sec. IIII Al and IIII Bl discuss 
the statistical properties of local concentration and tem- 
perature fluctuations, respectively. In Sec. IIII CI we show 
that there is a certain buoyant scale in turbulent thermal 
convection, above which buoyancy effects are important 
and temperature is active. Sec. IIII Dl compares the lev- 
els of intermittency and isotropy of the two scalars, both 
above and below the buoyant scale. Sec. IIII El studies the 
cross-correlation between the two scalars. We summary 
our findings and conclude in Sec. IV. 



II. EXPERIMENTAL SETUP AND 
PROCEDURES 

A. Convection cell 

The scalar measurements were conducted in a rectan- 
gular turbulent RBC cell, whose schematic drawing is 
shown in Fig. [T] The length, width and height of the 
cell are 25.4 x 7.5 x 25.4 cm^, respectively, and thus the 
large-scale circulation (LSCHs confined in the plane with 
the aspect ratio unity [H, [13, [III • The sidewall of the 
cell, indicated as E in the figure, is composed by four 
transparent Plexiglas plates of thickness 7 mm. The top 
(B) and bottom (F) plates, whose surfaces are electro- 
plated with Nickel to prevent the oxidation by water, are 
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FIG. 1: Schematic diagram of the rectangular cell, (a) Front 
view of the rectangular cell and (b) cross-section view of the 
bottom plate. A: the chamber cover, B: the top plate, CI, C2, 
C3, C4: nozzles, D: thermistors, E: the Plexiglas sidewall, F: 
the bottom plate, G: the cover for the bottom plate, Ifl: slit 
for dye injection, H2: nozzle for dilution water injection, H3: 
nozzle for inserting the small thermistor, H4: nozzle for fluid 
out-flowing, I: heater. 



made of copper with a thickness of 1.5 cm for its good 
thermal conductivity (400 W/mK). Rubber 0-rings are 
placed between the copper plates and the sidewall plates 
to avoid fluid leakage. A water chamber is constructed 
with the upper surface of the top plate and an attached 
stainless cover (A). There are four nozzles (CI, C2, C3 
and C4), located at the end of the long edges, through 
which a refrigerated circulator is connected to the cham- 
ber. To keep the temperature of the upper plate uniform, 
water is pumped into the chamber through two nozzles 
CI and C3, circulates within the chamber, and then goes 
out from the other two nozzles C2 and C4, respectively. 
A rectangular film heater (I) is sandwiched between the 
lower copper plate (F) and a stainless steel plate (G). 
There is a slit (HI), through which the fiuorescence dye 
solution is injected into the cell, parallel to the short edge 
and 1 cm from the sidewall on the lower copper plate. 
To minimize the perturbation of the bottom plate and 
the lower boundary layer, the area of the slit for injec- 
tion Ssiit is very small, which is 1 mm in width and 2 
cm in length, so the slit takes up only 0.1% of the bot- 
tom plate's total area. Diluted water is injected from 



3 



dilution water 





Convection Cell 
Injected dye 



Pre Amp 




DSA 






FIG. 2; (Color online). Schematic diagram of the experimen- 
tal setup for simultaneous measurements of temperature and 
concentration. 



the nozzle H2 to keep the background concentration of 
the cell unchanged and the nozzle H4 is used for water 
out-flowing. 



B. Temperature measurements 



through a stainless steel capillary tubing (U143, Fisher 
Scientific) with outer diameter of 1.5 mm and inner di- 
ameter of 1 mm and then attached to a translational 
stage with an accuracy to 0.001 mm. This arrangement 
allows the precise adjustment of the separation between 
the temperature probe and the laser beam, which is used 
for the local concentration measurement and will be in- 
troduced in Sec. Ill C II A Wheatstone bridge with AC 
source and a lock-in amplifier (Model SR830 DSP, Stan- 
ford Research Systems) are used to measure the temper- 
ature fiuctuations. A Dynamic Signal Analyzer (DSA, 
HP 35670A) with a 24-bit dynamic range and four input 
channels is used to digitize and record the signals and 
then transmit them to computer. 



C. LIF technique 

There are many experimental techniques for measuring 
a passive scalar field in turbulent flows, among which the 
LIF (laser-induced fluorescence) technique seems to pro- 
vide the best non-intrusive method for the high-spatial- 
resolution measurement of concentration. The use of this 
technique in fluid mechanics have been well documented 
[Tg, 20, 21, 22, 23, 24, 25]. In almost of aU these studies 
were carried out in open systems, such as turbulent flows 
in a pipe or in a channel. Whereas, LIF experiments in- 
volving longtime measurements in closed systems, such 
as the turbulent RBC, have not been reported. Longtime 
measurements in a close system will result in an increase 
of mean concentration of the system and hence leads re- 
ductions in the dynamical range of the fluorescence signal 
and also in the signal-to-noise ratio. In this section we 
present the details of the apparatus for our LIF measure- 
ments and how we overcome such difficulties, and show 
the validation of this technique in our case. 



There are two different types of thermistors used in 
our experiments for the temperature measurements. Five 
thermistors of the first type (Model 44031, Omega En- 
gineering Inc.), indicated as D in Fig. [U are imbedded 
inside the two plates to monitor their respective temper- 
atures. Three thermistors are in the top plate and the 
other two are in the bottom one. A 6.5-digit multime- 
ter (Model 2000, Keithley Instruments Inc.) is used to 
record the resistance of the thermistors with a sampling 
frequency of ~ 1 Hz. The measured relative temperature 
difference between two thermistors in the same plate is 
found to be less than 2% of that across the convection cell 
for both plates in our experimental conditions, indicat- 
ing that the temperature is uniform across the two hor- 
izontal plates. The second type of thermistors (AB6E3- 
B10KA103J, Thermometries) are used in the local tem- 
perature measurement inside the convecting fiuid. They 
have a sensing head of 250 /im in diameter and a ther- 
mal time constant of 10 ms in water, which is fast enough 
to detect the temperature fluctuations in thermal turbu- 
lence. As shown in Fig. [21 the thermistor is threaded 



1. Local concentration measurement 

Figure [5] shows the schematic diagram of the experi- 
mental setup for local concentration measurement. The 
solution of the fluorescence dye of Rhodamine 6G is in- 
jected from the bottom plate of the convection cell, and 
then a laser light source is used to illuminate the flow. 
From the absorption spectrum of dye solution with con- 
centration of 1.33 X 10^^ M measured by a monochro- 
mator (Model 300i, Acton Research Corporation) and 
an ICCD (Model TE/CCD-1024-E/1, Princeton Instru- 
ments Inc.), one sees that there is an absorbtion peak of 
fluorescence dye around the wavelength of 520 nm. Thus 
the green light with the wavelength of 514.5 nm from 
an argon ion laser (Coherent Innova 70) is chosen as the 
exciting light. For low concentration, the resulting fluo- 
rescence intensity at one point is proportional to the lo- 
cal concentration. We use a photomultiplier tube (PMT, 
R2368, Hamamatsu) to quantify the intensity of the fluo- 
rescence. From the fluorescence spectrum of dye solution 
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measured by the same equipments as for the absorbtion 
spectrum, there is an emission peak around the wave- 
length of 560 nm, so a filter with a passband from 550 nm 
to 570 nm is chosen and is placed in front of the PMT to 
block out the incident green light. To improve the signal- 
to-noise ratio, the experiments were conducted in a dark 
environment to minimize noise due to ambient light and a 
low-noise preamplifier (Model SR560, Stanford Research 
Systems) is used to filter out high frequency noise and 
amplify the signal obtained from PMT. The four-channel 
DSA, which is also used for recording the temperature 
fluctuations, digitizes and records the fluorescence sig- 
nals and then transmits them to computer. 

Since the turbulent RBC is a closed system, a contin- 
uous injection of dye will increase the background con- 
centration level in the convection cell and thus reduce 
the dynamical range of the fluorescence signal. To solve 
this problem, we injected dilution water into the cell to 
keep the background concentration constant. As shown 
in Fig. [21 both the dye solution and dilution water are 
supplied to the cell from the raised tanks and by the pres- 
sure difference between the tanks and the cell and hence 
the flow rate of injected fluids can be varied by changing 
the height of the tanks. The out-flowing fluid is captured 
in a third tank (not shown in the figure) from the nozzle 
H4 in Fig. 1. To minimize the perturbation of the flow 
field in the convection cell, the speed of injected fluids is 
approximately the same as that of the large-scale circula- 
tion (LSC) of the convective flow, which is about 1 cm/s 
in our case. In the present case, the volume of injected 
fluids per unit time is very small (about 0.2 mL/s for the 
dye solution and about 0.4 mL/s for the dilution water, 
in contrast to the convecting fluid volume of 4.84 L). 



2. Calibration for LIF 

Previous works have pointed out that there are some 
effects which can influence LIF measurement and hence 
chaiige the measured magnitude of the dye concentration 
[H, [2fl, m m, (Hi. The effects which we must consider 
in our experiments are as follows: 

PH and temperature efTects: The measured fluo- 
rescence intensity If may vary with PH and temperature 
[20| . However, in our experiment, the PH of the fluid 
in the convection cell is almost constant (^ 5.5) and 
the temperature fluctuations at the measuring position 
is about 2°C. In this range, PH and temperature effects 
can be neglected. 

Absorption along the illumination path: For 
our optical configuration, the incident laser must pass 
through the background fluid which itself contains dye 
in it before illuminating the measuring position. This 
causes the attenuation of the illumination source along 
its path due to absorption and then a non-linear relation- 
ship between fluorescence intensity and dye concentra- 
tion (23, ll^l . However for sufficiently low concentration, 
absorption effects along the incident laser beam path can 
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FIG. 3: (a) A sample time series of if{t) measured at the flow 
ceU center. Data taken at C = 2.0 x 10"^ M, le = 200 mW 
and V = 0.14 cm/s. (b) // (= {if{t))) vs v at h = 200 mW 
for, from bottom to top, C = 0.5 x 10"'' M, 1.5 x 10"'^ M, 
2.5 X 10"'^ M, and 4.0 x 10"^ M. 



be neglected and the fluorescence intensity at the mea- 
suring position can be expressed linearly as 



(1) 



where /3 is a proportionality constant, is the intensity 
of the induced laser and C is the dye concentration. 

Photobleaching: If the exciting laser intensity is 
high enough or the velocity v of the fluid is zero or too 
small for the LIF measurement, some dye molecules, af- 
ter excited, don't absorb or re-emit light and the fluo- 
rescence intensity decays with time, due to photodecom- 
position or collision quenching, etc. This phenomenon 
is called photobleaching and have been well studied in 
many experiments ^19,,^, ^4, 25]. 

To determine the proper ranges in which our experi- 
ments can be carried out properly, we calibrated the fluo- 
rescence intensity // under controlled conditions (e.g. C, 
le and v). The convection cell in Fig. [5] was replaced by 
a flow cell and a rectangular water tank, which contains 
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FIG. 4: (a) 7/ vs C at 7^ = 200 mW for three diflferent 
velocities. For clarity, the v = 0.28 cm/s and 0.57 cm/s data 
have been shifted up by 0.3 and 0.6, respectively, (b) 7/ vs 7e 
at u = 0.28 cm/s for two different concentrations. The solid 
lines show the linear fits to the corresponding data. 



dye solution with known concentration, is used as a reser- 
voir. Water was pumped from the tank to the flow cell 
and then returning back to the tank by a gear pump. The 
fluorescence intensity was measured at the center point 
of the flow cell. Assuming a laminar velocity proflle, we 
can obtain the mean velocity at the measuring position 
from the flow rate of the pump and the cross section area 
of the flow cell. 

A sample time-series of the measured instantaneous 
fluorescence intensity i f (t) is shown in Fig. EJa) . One 
sees no attenuation of if{t) in this 60-second segment. 
In Fig. ^h), we show // {— {if{t)), where (• ■ • ) denotes 
a time average) measured at different dye concentrations 
as functions of the flow velocity. We see no appreciable 
dependence of // on the fluid velocity v. In our scalar ex- 
periment, the mean velocity at the measuring position is 
about 1.8 cm/s at Ra = 1.2 x 10^° (see Sec. HiTeI . So in 
this range of fluid velocity photobleaching effects can be 
neglected. In Fig. Sfa) we show measured fluorescence 



intensity, at several flow velocities, vs the concentration 
C. It is seen that, over the range measured (<5.0xl0^^ 
M), // is linear with C. In Fig. H^b) // was measured 
with increasing incident laser intensity /g, which was 
measured by a power meter (model 1825-C, Newport). 
For le < 500 mW, // is found to be proportional to le- 
Based on the above measurements, we chose the concen- 
tration of the injected dye solution Cq to be 4.0x10"^ 
M, the background concentration Cback = 1.33 x 10"'' 
M and the incident laser intensity le — 200 mW in our 
scalar experiment. In these parameter ranges Eq. ([1]) is 
valid and the fluorescence intensity If is proportional to 
the local concentration C. Thus the local concentration 
can be obtained directly from the measured fluorescence 
intensity. 



D. Experimental conditions and parameters 

In our experiments, water is used as the convecting 
fluid and Rhodamine 6G as the fluorescence dye. During 
the experiments, the constant power is supplied to the 
bottom plate of the convection cell, so that it is under a 
constant-flux boundary condition, but at steady state its 
temperature remains effectively constant; the top plate's 
temperature is regulated so that it is under constant- 
temperature boundary condition. The control parame- 
ters are the Rayleigh number Ra, the Prandtl number Pr 
and the Schmidt number Sc, which are deflned as. 



Ra - 



a 



gATH^ 



, Pr 



, and Sc 



V 



(2) 



respectively. Here a, kq, and kt are the coefficients 
of thermal expansion, kinematic viscosity, molecular dif- 
fusivity of the fluorescence dye and thermal diffusivity of 
the working fluid, respectively, g is the acceleration due 
to gravity, AT is the temperature difference between the 
bottom and the top plates of the convection cell, and H 
is the height of the cell. Another parameter used to char- 
acterize the rate of advection of the scalar by a flow to 
its rate of diffusion is Peclet number, which is deflned as, 



uH , uH 
Pct — and Pec — 

KT KC 



(3) 



for temperature and concentration, respectively. Here u 
is the mean velocity at the measuring position. Table I 
summarizes the typical experimental parameters of the 
present investigation. From the table, one sees that, al- 
though Pec is two orders larger than Pct, both Pec and 
Pex are much larger than 1, suggesting that the mixing of 
the temperature and concentration flelds are both dom- 
inated by the turbulent velocity field, rather than their 
diffusion rates. Table I also shows that the temperature 
difference across the cell AT is 27.85°C. Previous studies 
[2^ [27} have pointed out that to strictly conform to the 
Boussinesq condition (in water) AT should be limited to 
< 15 °C. Thus, there may exist some Non-Boussinesq 
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TABLE I: Typical experimental parameters, u is determined from the cross-correlation between temperature and concentration 
fluctuations in Sec. IIIIEI ts and to (— 4:H/u) are the Bolgiano time scale and the turnover time of LSC, respectively. 



AT 


H 


Ra 


Pr 


So 


Per 


Pec 


Nu 


(°C) 


(cm) 














27.85 


25.4 


1.2 X 10'" 


5.3 


2100 


3.16 X 10** 


1.25 X 10' 


140 


Tbuik 




Kt 


KG 


u 


tB 


to 




(°C) 




(m^/s) 




(cm/s) 


(sec) 


(sec) 




30.97 


7.85 X 10-' 


1.48 X 10- ' 


3.73 X 10-"' 


1.84 


1.3 


55.2 





effect in our current system. Non-Boussinesq effect in 
water manifests as an increase in tlie mean bulk temper- 
ature as compared to the average of the top and bottom 
temperatures, and it also shghtly reduces the heat flux, 
i.e. Nusselt number [2^. However, its effect on the tem- 
perature and concentration fluctuations are unknown. 

In addition to the injection from the bottom plate, 
the dye solution was also injected into the central region 
of the convection cell. Due to the absence of the mean 
flow in the cell center, however, the dye packets meander 
around after injection and we could not acquire sufficient 
amount of mixing data (signals) for statical analysis. The 
measurements were also carried out near the sidewall of 
a cylindrical cell. Comparable features and behaviors 
of the passive scalar and of its relations with the active 
scalar were obtained. However, due to azimuthal mean- 
dering of the LSC in the cylindrical cell 29, 30, iii], the 
injected dye solution wobbles with the LSC and cannot 
always pass through the laser beam, resulting in a poor 
signal-to-noise ratio. The statistics are thus worse for the 
cylindrical cell. Therefore, all results presented and dis- 
cussed in this paper are from measurements made near 
the sidewall of the rectangular cell with dye injected from 
the bottom plate. 

For the present Ra and Pr, the viscous boundary layer 
thickness near the sidewall is about 2 mm [s^l • Thus the 
measuring position is chosen at 1 cm from the sidewall 
and at the mid-height of the convection cell, which is 
outside of the sidewall viscous boundary layer and inside 
the LSC. The cell is tilted with a small angle of about 
1° near the side of the dye injection slit to lock the LSC 
and to make the dye solution move upwards with the 
LSC after it is injected into the cell. 

To reveal the statistical properties of passive and ac- 
tive scalars, longtime measurements of the local concen- 
tration and temperature fluctuations were carried out in- 
dependently. The measurement of concentration lasted 
27 hours with a sampling frequency of 64 Hz. While 
the time series of temperature fluctuations, which was 
obtained at the same measuring position and under the 
same conditions as concentration, consists of three parts, 
each of which lasted 20 hours and all were acquired at a 
sampling rate of 64 Hz. To study the relations between 
these two scalars, a simultaneous measurement of the lo- 
cal concentration and temperature fluctuations was also 
made. In this case, the temperature probe (thermistor) 
was placed above the concentration probe (laser beam). 
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FIG. 5: (Color online), (a) A 60-second time series of concen- 
tration fluctuations, (b) PDF of concentration fluctuations 
P{C). The dashed curve represents a Gaussian-function flt- 
ting to the background signal and the solid line is for reference. 



The separation £ between two probes can be adjusted by 
a translational stage. For each separation, we recorded 
two hours of data, except the measurement for £ = 0.375 
mm which lasted 6 hours, with the sampling frequency 
of 128 Hz for both concentration and temperature. 
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FIG. 6: Power spectra -Pa(/) of the fluctuating concentration 
C{t). Inset is the plot of log[Ps(/)] as a function of j"'*^ for 
the same data. The solid line is for reference. 



Figure [6] shows the power spectra Pg (/) of the mea- 
sured dye concentration C. The figure shows that the 
magnitude of Ps{f) spans about five decades from the 
large scale to the noise level. And we see no obvious 
scaling in the low- frequency range, which may be due 
to a low Reynolds number (about 5700 in our scalar ex- 
periments). At higher frequencies, the measured Ps{f) 
drops sharply and can be fitted by a stretched exponen- 
tial function 



P^(/) ~ e[^(^-^/-^'=)''], 5 = 0.49 ±0.05 



(4) 



where fc is a cutoff frequency. The inset of Fig. [S] shows a 
semi- logarithmic plot of the measured Ps{f) as a function 
of In this plot, the stretched exponential function 

appears as a straight line. The cut-off frequency of the 
concentration is around 10 ~ 20Hz, which approximately 
equals to that of the temperature (see Sec. IIII Bl) . 



III. RESULTS AND DISCUSSION 

A. The passive scalar measurement 

Figure [HKa) shows a typical time series of the mea- 
sured concentration fluctuations. When the injected dye 
packets passing through the concentration probe, a sig- 
nal with a pulse-like shape would be detected. From the 
figure, one sees that our measured C{t) consists of many 
such pulse-like signals of variable heights, corresponding 
to dye packets with different concentration levels. 

In homogeneous turbulent flows, probability density 
function (PDF) of the velocity fluctuations is Gaussian 
[ssj (It was also showed that the velocity fluctuations in 
turbulent RB system is Gaussian distributed [H, Issj). 
It is generally assumed that the PDF of a passive scalar 
is also Gaussian, since it is advected by the Gaussian- 
distributed turbulent velocity field. Earlier measure- 
ments about the scalar PDF in homogeneous turbulence 
p}6j showed that a Gaussian distribution was a satisfac- 
tory model for the scalar signal itself except that the tails 
of the PDF have not been examined. However, Jayesh & 
Warhaft found that for a mean temperature gradi- 
ent in grid turbulence the tails of the temperature PDF 
have an exponential shape (in their system temperature 
is passive). Our measured PDF P{C) is shown in Fig. 
[S^b). One sees that the PDF consists of two parts. The 
left part, which has a abscissa value below —0.25, is the 
distribution for the background signals. The distribution 
of this part is approximately Gaussian, and one can use a 
Gaussian function to fit the data of this part (the dashed 
curve in the figure). The right part is the distribution of 
the mixed fiuorescence dye. It can be seen that this part 
exhibits an exponential decay, which was also observed 
for passive scalars in other types of turbulent fiows, e.g. 
in a turbulent jet [37| . 



B. The active scalar measurement 

Figure EKa) shows how the temperature fluctuations 
change with time. The measured T{t) shows many inter- 
mittent upward sharp spikes of variable heights, which 
may be attributed to the upward-moving hot plumes. 
The corresponding PDF of the normalized temperature 
fluctuations is shown in Fig. [7fb). Unlike the expo- 
nential distributions observed in the center and near the 
sidewall of a cylindrical cell, the measured PDF near the 
sidewall in the rectangular cell shows a flat central re- 
gion and two Gaussian-like tails, which are fltted by sep- 
arate Gaussian functions (the solid and dashed curves). 
The Gaussian-like distribution of temperature fluctua- 
tions means that the temperature probe felt a lot of ther- 
mal coherent structures, i.e. thermal plumes, passing 
through it [s^]. The ratio of the long side to the short 
side of our rectangular convection cell is about 3.3, so 
that the LSC is largely locked in the vertical plane paral- 
lel to the long side. Hence most of thermal plumes would 
group together, move upwards with the LSC and pass 
through the measuring position due to the conflnement 
in the other direction. Thus, although our experiments 
with Ra = 1.2 x 10^° were carried out in the so-call "hard 
turbulence" regime [sj, the measured temperature fluc- 
tuations near the sidewall have two Gaussian-like tails 
rather than exponential ones. It is also not surprising 
that the right peak is more than ten times higher than 
the left one, since the LSC has been locked and at the 
measuring position there are more hot plumes moving up- 
wards than cold one moving downwards. Note also that 
there is only one group of cold plumes in the 60-second 
time series shown in Fig. [7l[a), as compared to the large 
number of hot plumes in the same plot. 

We now turn to the power spectra of temperature fluc- 
tuations. Based on dimensional arguments, it is tradi- 
tionally thought that there is simple scaling relation be- 
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FIG. 7: (Color online), (a) A 60-second time series of tem- 
perature fluctuations. The dashed line indicates the mean 
bulk temperature Ttuik- (b) Corresponding PDF of normal- 
ized temperature fluctuations P(T). The dashed curve and 
the solid curve are two Gaussian-function flttings to the left 
and right tails, respectively. 
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FIG. 8: Power spectra Ps{f) of the fluctuating temperature 
T{t). Inset is the plot of \og[Ps{f)] as a function of f '"^ for 
the same data. The solid line is for reference. 



tween Ps{k) and k in the low-frequency region, i.e. 



(5) 



where k is the wavenumber. Here, a would be 5/3 ac- 
cording to Kolmogorov's arguments (K41)[39l[4oL [ 4l | .l4^. 
and it would be 7/5 if the arguments of Bolgiano [43| and 
Obukhov [4J] (B059) are correct. However, from a direct 
multipoint measurement of the temperatrue field. Sun et 
al. [45] showed that a = 8/5 near the sidewall of a cylin- 
drical cell, which could be explained by a simple scaling 
analysis based on the coaction of buoyancy and inertial 
forces. 

Figure [H shows the measured temperature power spec- 
tra Ps{f) as a function of /. From the figure, one 
sees that the magnitude of Ps{f) spans more than nine 
decades from the large scale to the noise level. However, 
the scaling range in the low-frequency region is so lim- 
ited that one cannot tell the accurate value of a. This 
may be due to the asymmetry of the rectangular cell 
and the sidewall effect. Another possible reason for this 
limited scaling range is that the measured Ps{f) here is 
based on time or frequency domain. By using the Taylor 
frozen-turbulence hypothesis with the condition that the 
turbulent velocity fluctuations is much smaller than the 
mean flow velocity, one can relate the time domain re- 
sults to those of the space domain. However, it is known 
that the mean flow velocity is comparable to the rms ve- 
locity near the sidewall in turbulent RBC convection cell 
[4^. This would flatten the measured Ps{f) in the low- 
frequency region. At higher frequencies, the measured 
Ps{f ) drops sharply and can be fitted by a stretched ex- 
ponential function 



Psif) 



if/fcf 



6 = 0.70 ±0.05 



(6) 



where /c is a cutoff frequency. (We also used the data 
measured in a 81 x 20 x 81-cm'^ rectangular cell from 
Ref. [TtI to study the temperature power spectra near 
the sidewall and b = 0.73 ± 0.05 was obtained.) Previous 
measurements of temperature in smooth [l^l and rough 
[m cylindrical cells also showed a stretched exponential 
decay at higher frequencies. However, the values of b in 
those studies are both 0.55 ±0.05, which are smaller than 
our results. Power spectrum of the temperature fluctua- 
tions describe the cascade of entropy of the system and 
hence a larger value of b here means a more rapid cas- 
cade of the system's entropy compared with the cylindri- 
cal convection cell. This may be due to the asymmetry 
of the rectangular cell and the sidewall effect. The inset 
of Fig. [5] shows a semi-logarithmic plot of the measured 
Psif) as a function of Z^ ™. In this plot, the stretched 
exponential function appears as a straight line. From 
Fig. [51 one sees that the cut-off frequency of the temper- 
ature is around 10 ~ 15Hz, which is about the same as 
that of concentration as discussed in Sec. IIII Al 
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FIG. 9: (Color online). Cross-correlation coefBcient Rq'^{t) 
of the concentration and temperature increments. The two 
dashed lines show tB and to. The separation i between the 
temperature probe and the laser beam is 0.375 mm. 



C. The buoyant scale of convective turbulence 



In 1959, Bolgiano [43| and Obukhov [4j| proposed that 
there is a characteristic length scale in buoyancy-driven 
turbulence, now commonly referred to as the Bolgiano 
scale £b, above which buoyancy effects are important. 
It has been shown [49| that the Bolgiano scale can be 
estimated as 



5/4 



3/4. 



(7) 



where and ex are the kinematic and thermal dissipa- 
tion rates. Note that if the local dissipation rates e„(r) 
and CTir) are used in the above equation, then it gives a 
local Bolgiano scale. The global Bolgiano scale may be 
obtained by using the expressions of the globally aver- 
aged dissipation rates and in terms of Ra, Pr, and 
Nu which gives: 



(8) 



As all our measurements are made in the time domain, 
we use ts, which is the counterpart of £3 in the time 
domain, as a characteristic time scale in our system, i.e. 



tB = 



(9) 



(i?aPr)i/4' 

which can be easily evaluated from the measured values 
of Ra, Pr, Nu and the large-scale flow turnover time t^. 

To check whether this buoyant scale exists in turbu- 
lent RBC, we compute the cross-correlation coefficient 
Rq'^{t) between the increments of concentration and 
temperature, i.e., 

{drC{t)SrT{t + e/u)) 



r 



(10) 



where SrC(t) = C{t + T)-C{t), 5rT{t) = T{t + T)-T{t), 
and £/u describes the time delay between the concentra- 
tion and temperature signals due to the separation of the 
two probes. When the concentration and temperature in- 
crements are uncorrelated, we have -Rjr'^(T) = 0, while 
Rq^{t) = 1 is obtained for a linear relation between the 
two increments. 

Figure [9] shows the measured Rq'^{t) as a function of t 
for £ — 0.375 mm, which is the measured separation be- 
tween the laser beam and the thermistor tip. Two dashed 
lines in the figure show ts and to, which were determined 
from the measured Ra, Pr, and Nu (see Table I) using Eq. 
(9) and from AH/u respectively. One sees that Rq-^{t) 
increases with increasing t when t is small and saturates 
at a value of ~ 0.45 when t > ts- It has been shown 
previously [5Q[ that the vertical velocity increments cor- 
relate strongly with the temperature increments above 
the Bolgiano time scale tB, reflecting strong effects of 
buoyancy. Since the concentration fluctuations are dom- 
inated by the vertical velocity field near the sidewall, a 
large positive correlation between the concentration and 
temperature increments is expected and indeed observed 
here. This supports the idea that, above the Bolgiano 
time scale, buoyancy effects are significant and tempera- 
ture can be regarded as an active scalar. In contrast, for 
time scales below tB , buoyancy is not important and the 
temperature field behaves as a passive scalar. Thus, the 
concentration and temperature fluctuations are two sep- 
arate stochastic fields at scales below tB- In this case, the 
correlation between them should be very weak (as shown 
in Fig. [HI Rq^{t) decreases sharply when r < tB), even 
though they are advected by the same velocity field. 

The physical picture revealed by Fig. [5] is that a cer- 
tain buoyant scale exists in turbulent RBC, above which 
buoyancy effects are important. Thus, the temperature 
field can indeed be regarded as an active scalar above tB, 
and as a passive scalar below tB- 



D. Statistical properties of scalar fluctuations 



We have seen in Sec. IIII Al that the right tail of 
concentration PDF, which describes the distribution of 
the mixed fluorescence dye, has a decreasing exponen- 
tial shape. An exponential distribution near the tail 
means higher probability of rare events, i.e. intermit- 
tency, in comparison with the Gaussian distribution. The 
small-scale intermittency of concentration fluctuations 
can be characterized by the distributions of their in- 
crements over different scales. Figure [TUla) shows the 
measured PDFs of the concentration increments ST-C{t) 
for four different time lags r, with r = 0.05tB (close 
to the dissipative scale), 0.15*^, tB (Bolgiano scale), 
and lOtB to/A), from top to bottom. The PDFs 
have been normalized by their own standard deviations 
= (('^tC')^)^/^ and shifted vertically for clarity. The 
dashed line in the figure represents a Gaussian distribu- 
tion of variance 1 for reference. Three features are worthy 
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FIG. 10: (Color online). PDFs of (a) concentration and (b) 
temperature increments normalized by their standard devi- 
ations for four time intervals r/tB ~ 0.05, 0.15, 1, and 10, 
from top to bottom. The curves have been shifted vertically 
for clarity and the black dashed curve indicates a Gaussian 
distribution of variance 1. 



of note: (i) all PDFs are strongly non-Gaussian from the 
smallest time interval to lags of the order of the large 
scale; (ii) all PDFs are clearly not self-similar, indicating 
strong intermittency; (iii) there is a slight asymmetry of 
the PDFs' tails, especially for small time lags, reflect- 
ing the well-known property of small-scale persistence of 
anisotropy, i.e. ramp-chff structures [H, H, 0]. 

The distributions of the normalized temperature incre- 
ments S-rT/ar (ct = {{^rT)^y^^) for the same values of 
T as those for the concentration are also examined, and 
are plotted in Fig. [TPl b). The small-scale persistence of 
anisotropy can also be found for the active scalar incre- 
ments from their asymmetric distributions at small time 
lags. This asymmetry can be linked to the sharp fronts 
of coherent structures in turbulent thermal convection, 
i.e. thermal plumes [lB|. Another notable feature is that 
as time lags increase towards the large time scale the 
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FIG. 11: (Color online), (a) Skewness and (b) flatness of 
the concentration (red circles) and temperature (blue trian- 
gles) increments as functions of t/Ib- Inset of (b): the local 
slopes, X^i'^/^B) = dlog (Fx(r/ts))/dlog (r/ifl), as func- 
tions of r/ts- 



shapes of the PDFs evolve from approximately exponen- 
tial tails to nearly Gaussian-like distribution, reflecting 
strong intermittency as well. 

When comparing the PDFs of concentration and tem- 
perature increments, two features are remarkable: (i) 
The PDFs of temperature increments at small scales [the 
top two curves in Fig. llOf b)] are more asymmetric than 
those of concentration increments at the corresponding 
scales, suggesting that the mixing of temperature are 
more anisotropic than those of concentration, (ii) Qual- 
itatively, Fig. [TO] shows that the PDFs of temperature 
increments change from non-Gaussian in small scales to 
near Gaussian in large scales, whereas, those of concen- 
tration increments remain non-Gaussian for all scales in- 
vestigated. In this sense, it suggests that the shapes of 
the concentration increments' PDFs are less dissimilar 
among themselves than those of the temperature incre- 
ments' PDFs as time lags t increase towards the large 
time scale, and hence the passive scalar is less intermit- 
tent than the active one [51| , which can also be quantified 
by the flatness discussed below. 
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To quantitatively characterize the anisotropic effects of 
the small-scale mixing of the two scalars, we study the 
skewness of their increments, i.e., 



Sxir) = {{SrXr)/{{SrXrf/\ 



(11) 



where SrX{t) = X{t + t) - X{t) with X = C for con- 
centration and X = T for temperature. The skewness is 
a global measure of asymmetry. By definition, the skew- 
ness of a symmetric-distributed quantity is 0. Figure 
ITlT a) shows the measured skewness for the concentra- 
tion and temperature increments. One sees that both 
Sc{t) and St{t) decrease with time lags for small r. 
Sc{t) reaches its minimum value around r ~ O.T^b and 
then increases a little, whereas, 5't(t) decreases almost 
monotonically to the zero value. For t < ts, both tem- 
perature and concentration behave as passive scalars but 
one finds \St{t)\ > \Sc{t)\, suggesting that the mixing 
of temperature are more anisotropic than those of con- 
centration. This is because the mixing of a passive scalar 
field becomes more isotropic with decreasing diffusivity 
[5^ and from Table I we see that the thermal diffusivity 
kt is about two orders of magnitude larger than molec- 
ular diffusivity kc of the dye. For t > Ib, both Sc{t) 
and 5't(t) fluctuate around 0, implying that the mixing 
of both scalar are isotropic above the buoyant scale. This 
is strange since buoyancy effects act on fluids along the 
vertical direction but do not induce anisotropy for both 
temperature and concentration. 

We now examine more quantitatively the intermit- 
tency levels of the two scalar flelds. In general, the 
deviations of scaling exponents of structure functions, 
Rp (t) = {\6tXY') ^ T^p , in the inertial range from those 
expected by the arguments of Kolmogorov [H, \^ or Bol- 
giano [ij and Obukhov [l^l can be used as a measure 
for the degree of intermittency. Unfortunately, above the 
Bolgiano time scale is , no discernible scaling for concen- 
tration is observed no matter whether structure functions 
are plotted directly against time lags in a log-log plot 
or indirectly using an extended self-similarity method 
[53L 113] ■ Here, we use the flatness of the temperature 
and concentration increments, defined as, 

Fxir) = {{5,Xf)/{{6^Xff = Rf /{R^f - 

(12) 

to characterize their levels of intermittency, where — 
— 2£^2 ■ The fiatness is a characteristic measure of 
whether the data are peaked or flat relative to a Gaus- 
sian distribution and hence can be used as a relatively 
sensitive measure of intermittency. By definition, the 
fiatness of a Gaussian-distributed quantity is 3. In addi- 
tion, the local slope of the flatness for the temperature 
and concentration increments, defined as 



X^(t) =dlog(Fx(T))/dlogr, 



(13) 



may also be used as a measure of self-similarity, since 
X"^(r) = within a certain range of scales for a quan- 
tity whose increments are scale-free within that range. 



Therefore, larger deviations of (''") from reflects less 
self-similarity, i.e. more intermittent. Figure fTU b) shows 
the measured Fc{t) and Ft{t) plotted as functions of 
time lags r. One sees that both quantities decrease with 
increasing r for small r and cross at r ~ O.Sis. For 
r > lots, Ft{t) levels off to a value of about 3.5, which 
is very close to that of a Gaussian-distributed quantity, 
while Fc{t) levels off at a value of 5 for r > O.S^b- The 
local slopes of Fc{t) and Ft{t) are plotted as functions 
of T in the inset of Fig. [TlT b). Above ts, although 
one finds that the value of Ft{t) is smaller than that 
of Fc{t), temperature possesses a higher level of inter- 
mittency. This is because x'^i''') — 0, suggesting an ap- 
proximate self-similar property, and |x^(''')l > \x'^{T)\y 
i.e. the deviations of x^(''') from are larger than those 
of X*^ (''')■ Recall that within such regime buoyancy is 
important and the temperature field is an active scalar. 
Therefore, our results presented here suggest that the 
active scalar has a higher level of intermittency than 
the passive scalar, at least in turbulent thermal con- 
vection. Below ts, the situation is more complicated. 
For O.lSts < T < ts, \x^{t)\ > \x^{t)\ is obtained, 
again suggesting that temperature possesses a higher 
level of intermittency within such regime. Whereas, for 
T < O.lSis, we have |x"^(''')l < \x'"{''')\^ suggesting that 
concentration is more intermittent for such regime. Note 
that this regime is close to the dissipative scale, but we 
don't understand the mechanism that caused such result. 



E. Cross-correlation between the concentration 
and temperature fluctuations 

Figure [TWa) shows the measured auto-correlation 
functions, 



{6X{t)SX{t + T)) 



(14) 



for temperature {X = T) and concentration {X = C), 
respectively, and their cross-correlation function 



9c[T) 



{5C{t)5T{t + r)) 



(15) 



rms J- rms 



as functions of time lag r, where 5X{t) = X{t) — {X{t)). 
One sees that all three functions have a shape similar to 
each other, but the full width at half maximum (FWHM) 
of their peaks are different. The FWHMs of the temper- 
ature and concentration fluctuations are Wt — 3.7s and 
Wc — 0.5s, respectively, and the FWHM of their cross- 
correlation is Wtc — 2.2s, which is about the average of 
Wt and Wc- One half of the FWHM of the correlation 
functions, Wt/2 (iyc/2), provides information about the 
transient time for the temperature (concentration) fluc- 
tuations of a certain size passing through the temper- 
ature (concentration) probe with a certain speed, and 
thus one can obtain the coherent length of the temper- 
ature (concentration) fluctuations from the FWHM and 
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FIG. 12: (Color online), (a) 
Autocorrelation functions 
glij) (red line), g^ij) (blue 
line) and cross-correlation 
function (?c(t) (black line). The 
separation between the laser 
beam and the temperature 
probe is ^ = 0.375 mm. (b) 
Cross-correlation functions 
(?c(t) for five different separa- 
tions i. (c) The separation I 
vs the position of the measured 
correlation peak tq. The solid 
line shows the fitting linear 
function i = mtq with u = 18.4 
mm/s. (d) The measured 
correlation amplitude gaif) as 
a function of the separation 
t. The solid curve shows the 
fitting exponential function 

go(^) = 0.429e~^/^, with 
C = 27.3 mm. 



the advected speed. Since both scalars are transported 
by the same velocity field and the relation Wt > Wc 
is obtained, we conclude that the coherent length of the 
active scalar, i.e. the size of thermal plumes, is larger 
than that of the passive scalar, i.e. the width of the cliff 
structures of the concentration field [H, H, , which 

is probably due to the larger value of the thermal diffu- 
sivity. 

Figure [T2r b) shows the measured concentration- 
temperature cross-correlation functions for five different 
sparation £. It is found that the measured g{T) all have 
a single peak go as concentration and temperature are 
both advected by the same velocity field. However, the 
peak position tq for different separation £ are not the 
same. The bigger £ is, the larger tq becomes. How the 
measured tq varies with the separation £ is plotted in Fig. 
[T2y c). It can be clearly seen that £ increases with increas- 
ing To, and the increasing manner may be described by 
a simple linear function 

£ = utq, u = 18.4 mm/s (16) 

The solid line in Fig. W2\ c] shows the fitting linear equa- 
tion (fT6|l . The slope u = 18.4 mm/s, describing the 
mean velocity with which the temperature and concen- 
tration fiuctuations pass through the measuring points. 
Fig. fTSFd) shows the relationship between the obtained 
(7o and the separation £. It is clearly shown that go(^) de- 
creases with increasing £, and the decay may be described 
by a simple exponential function 

goW = 0.429e~^/^, ,^ = 27.3 mm (17) 

The solid curve in Fig. [T2t^d) shows the fitting equa- 
tion (fT7|) . A similar correlation length was also obtained 
in velocity-temperature cross-correlation functions in a 
cylindrical cell [56l| . 



IV. SUMMARY AND CONCLUSION 



To summarize, we have carried out a comparative 
study on local mixing of active and passive scalars in 
turbulent thermal convection. The statistical properties 
of local temperature are experimentally compared with 
those of local concentration. It is found that, although 
transported by the same turbulent flow field, the fluc- 
tuations of the two scalars are distributed differently: a 
decreasing exponential tail is found for the fluctuating 
concentration, while Gaussian-like tails for the fluctuat- 
ing temperature. It is also shown that the temperature 
increments exhibit a positive and large correlation with 
the concentration increments above the Bolgiano time 
scale tB\ whereas this correlation decays with decreas- 
ing time lags r for t < <b. This reflects the fact that 
there is a certain buoyant scale in turbulent RBC, above 
which the buoyancy effects are important. Above is, 
temperature is active and the mixing of both scalars are 
isotropic even though buoyancy acts on the fluid in the 
vertical direction. While temperature fluctuations are 
found to be more intermittent than concentration fluc- 
tuations, suggesting that the active scalar possesses a 
higher level of intermittency in turbulent thermal con- 
vection. Below is, temperature is passive but its mixing 
are more anisotropic than those of concentration, reflect- 
ing the fact that its diffusivity (kt) is larger than that 
of concentration (kc)- It is further found, from the si- 
multaneous measurements of temperature and concen- 
tration, that the two scalars have similar autocorrelation 
functions and there is a strong and positive correlation 
between them. Future investigations will be focused on 
statistical and geometric features of the passive scalar 
field for varying Ra. 
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